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METHOD AND APPARATUS FOR CALCULATING 
BIT LOG-LIKELIHOOD RATIOS FOR QAM SIGNALS 



FIELD OF THE INVENTION 
The present invention relates generally to communications systems, and 
more particularly to methods for demodulating and decoding encoded data bits 
transmitting in a code division multiple access (CDMA) system. 



BACKGROUND OF THE INVENTION 
Various communication systems are known in the art. In many systems an 
information signal is modulated onto a carrier signal and transmitted from a first 
location to a second location. At the second location, the information signal is 
demodulated and recovered. 

Typically, the communication path used by such systems has various 
limitations, such as bandwidth. As a result, there are upper practical limitations that 
restrict the quantity of information that can be supported by the communication path 
over a given period of time. Various modulation techniques have been proposed that 
effectively increase the information handling capacity of the communication path as 
measured against other modulation techniques. 

One modulation technique is M-ary quadrature amplitude modulation 
(QAM). QAM provides a constellation of an M number of modulation values 
(distinguished from one another by each having a different combination of phase and 
amplitude) wherein each constellation point represents a plurality of information 
bits. In order to achieve greater spectral efficiency in CDMA systems, for example, 
higher order modulations are used. In general, the number of bits that are 
represented by each symbol in an M-ary QAM system is equal to log 2 M. For 
example, 16 QAM uses 16 distinct constellation points (four points in each quadrant 
of a complex Cartesian plane having I and Q axes) that each represent log 2 16 or four 
information bits. Another example is 64-QAM wherein 64 constellation points each 



represent a combination of log 2 64 or six information bits. 

FIG. 1 illustrates a constellation for a 64-QAM communication system that 
is a map of 64 points each representing a combination of six bits on the complex 
plane defined by a horizontal axis representing the real portions (i.e., I) and a 
vertical axis representing imaginary portions (i.e., Q) of a complex number. 
Transmitted QAM information symbols on a communication channel (and the pilot 
and sync symbols as well) are discrete packets of a carrier signal modulated to 
convey information using both the amplitude and phase-angle displacement of the 
carrier from some reference. QAM information symbols are represented on the 
constellation of FIG. 1 as complex quantities represented as vectors having both 
magnitude (represented as length or distance from origin) and phase angles $ 
(angles measured with respect to one of the axes). In a 64-QAM system, having 64 
different magnitude and phase angle combinations that correspond to 64 different 
possible bit patterns of six binary digits (which bits are from a serial stream of bits 
from an information source), each of the 64 points on the constellation is identified 
as representing one combination of six bits. 

A vector 100 (expressed in rectangular coordinates as 51 +5Q and having a 
length = (5 2 +5 2 )/2 and a phase angle 102 equal to the arc tangent of 5/5 or forty- 
five degrees with respect to the real axis), points to the point 104 in the 
constellation. This point is shown in FIG. 1 as representing the series of six binary 
digits (i.e., 000000). Another QAM symbol represented by point 106 (-11 - 1Q) is 
illustrated in this constellation and represents another series of six binary digits 
(i.e., 110011). 

For any modulation scheme, the demodulator preferably should provide the 
decoder with the log-likelihood ratio (LLR) of the received encoded bits. As shown 
in FIG. 2, a receiver 200 in a CDMA system includes a log-likelihood ratio 
calculator 202 that receives demodulated symbols from a demodulator 204 and 
calculates a LLR for each bit of the demodulated symbols. The decoder 206 
receives the LLR's calculated by the log-likelihood ratio calculator 202 and decodes 
the individual bits based on the received LLR's. 

When two-dimensional modulations (i.e., having I and Q components) yield 



spectral efficiencies greater than 2bps/Hz, each encoded bit cannot be mapped to an 
orthogonal symbol component. As a result, bit LLR's must be calculated from 
Euclidean distances between the demodulated symbol and the constellation symbols. 
The LLR of each actual j ,h bit associated with a symbol y received at time k 
(referred to as y t ) is determined by the following relationship: 



where u k j is the hypothesis j" 1 bit value of the symbol u at time k and P 
represents the probability that the f bit of u k is a binary value of 0 or 1 given that y k 
is the received symbol. This equation reduces to the expression: 



assuming that all transmitted symbols are equiprobable. The probability P is 
represented by the following relationship assuming a perfect channel correction, 
such that the real and imaginary components are independent: 




(1) 




(2) 



P(y k | «* ) = piy\ » y'k \ < > < ) = pM I u l ^O* I < ) 



(3) 



Assuming additive Gaussian noise, the above equation reduces to 
following expression: 
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where D\ is the squared Euclidean distance between the received symbol 
and the hypothesis symbol and o 2 is the variance of the Gaussian noise. 
Accordingly, the squared Euclidean distance D 2 k is calculated by summing the 
squares of the differences between the received and hypothesis symbol of both the 
real (r) and imaginary (i) components as represented by the following relationship: 



A 2 =k-»,r=w-<) 2+ <^-^ )2 (5> 

Substituting equation (4) into equation (2) allows calculation of the LLR 
based on the squared Euclidean distances and yields the following expression: 

-D\ -D 1 

LLR(u kJ ) = log Xe 2ff2 -log Yu e2 ° 2 (6) 

u k :u kJ =0 u k :u kJ = \ 

This expression may be approximated by simply taking the difference 
between the minimum squared distances between the received symbol and the 
hypothesis symbol having respective bit values of 1 and 0. Thus, equation (6) may 
be reduced to the following expression to approximate the LLR: 



LLR(u kJ )= minj^ 2 ] - min^D]} 



(7) 



where the LLR of a bit u kJ in a symbol is found by computing the difference 
of the minimum distance between the received symbol and a QAM constellation 
point where that particular bit equals one and the minimum distance between the 
25 received symbol and a QAM constellation point where that particular bit equals 



zero. In M-ary QAM systems the M number of symbols is equal to 2 m with m 
being the number of bits per symbol. Hence, a m number of LLR calculations are 
required per modulation symbol (i.e., finding the LLR for each of the m number of 
bits in the symbol). In order to determine the m number of LLR calculations, a M 
number of squared Euclidian distances (D 2 ) must first be determined between the 
received QAM symbol and each of the QAM symbols in the QAM constellation of 
symbols. Given the above mathematical operations, the calculation of LLR values 
for each received modulation symbol require a M number of squared Euclidean 
distance calculations (i.e., D 2 calculations, each of which involve 2 multiplication 
operations and an addition operation), a M number of comparison/selections per bit 
to determine the minimum squared Euclidean to the particular bit values of one and 
zero for all of the m bits in the symbol (i.e. , the number of min D 2 values where u 
= 1 (M/2 comparison/selections and the number of min D 2 values where u = 0 
(M/2 comparison/selections), the total then multiplied by the m number of total bits 
in the symbol and a m number of addition operations for determining the difference 
between the values as shown in Equation (7) above. Accordingly, the total number 
of operations required per received QAM modulation symbol is 2M multiplication 
operations plus m«M comparison/selection operation plus a M + m number of 
addition operations. 

For purposes of illustrating reduced frame processing time, it is assumed that 
a clock rate for a CDMA system is 32 x 1.2288 MHz, or 39.3216 MHz, although 
the teachings of the present application are advantageous at any number of different 
clock rates. A multiplication operation typically requires 2 clock cycles to 
complete, addition operations 1 clock cycle and comparison/selection operations 2 
clock cycles. Thus, for example, the time required to process a frame with 64- 
QAM modulation containing 5376 symbols with 6 bits per symbol, would be 
approximately 150 ms. This time period becomes problematic because it is much 
greater than the 5 to 20ms frame turn around time required for high speed data 
channels. Techniques such as parallelization and pipelining may be used to reduce 
calculation time, but are nonetheless insufficient to reduce the receiver throughput 
rate down to the order of a frame turn around time in a high speed data channel. 



Accordingly, there is a need for reduction of computation time for LLR's to 
the order of a frame time in a receiver that bit decodes QAM modulated symbols, 
especially for higher orders such as 64-QAM. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a Cartesian graph showing symbols in a 64-QAM system; 

FIG. 2 illustrates a block diagram of a receiver in a communications system; 

FIG. 3 illustrates a Cartesian graph showing Karnaugh mapped, Gray coded 
constellation symbols in a 64-QAM communication system according to the 
teachings of the present invention; and 

FIG. 4 illustrates a schematic of a logic circuit employed to calculate LLR 
values for each I component bits in a 64-QAM communication system according to 
the teachings of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The above noted problem is solved by the method and apparatus of the 
present disclosure. Specifically, mapping the QAM signal with a specific mapping 
scheme allows reduced computation time for log-likelihood ratios for each of the 
bits in a QAM symbol through the elimination of calculation steps based on 
characteristics of the predetermined specific mapping scheme. Preferably, the 
signal is Karnaugh mapped and Gray coded, which introduces symmetry into the 
constellation such that calculation steps may be eliminated based on the symmetrical 
characteristics of the QAM constellation. 

FIG. 3 illustrates an exemplary 64-QAM square constellation with each of 
the 64 points corresponding to a particular bit sequence of the 6 bits. Each of the 6 
bit sequences is of the form i,q 1 i 2 q 2 i 3 q 3 labeled from most significant bit to least 
significant bit. Of significance is that each of the "i" bits are independent of the 
ordinate or Q axis and each of the "q" bits are independent of the abscissa or I axis. 
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For example, ,he -r bits change along ,he -r or abscissa ax,s of ,he Cartesian 
gra ph in FIG. 3 (i.e., along .he rows), bu, remain .he same for al. cons.ella.ion 
po,n.s in ,he »Q" or ord.na.e axis direc.ion (i.e., .hey change from column .o 
column in .he abscissa direction). For example, .he bi. i, in .he second column to 
t he righ, of .he Q axis are a„ zeros, whereas all of .he i 2 bi.s for .he symbols ,n .he 
third column from .he righ. of .he Q axis are a., ones. To further illustrate th,s 
pr„ P er,y of .he Karnaugh mapped QAM signal, .he lines shown a, .he ,o P of .he 
graph and ,he left side of .he graph indica.e subse, of constellation poin, where.n 
th e bi. is a one value whereas .he absence of a line indica.es a subse. of pom, 
wherein .he corresponding bi. is a zero value. For example, .he subse, of po,n.s 
where bi. i, is equal .o a value of one al, lie in .he left half of .he Cartesian plane 

plane. 

in ,he following discussion a TRUE region of the constellation for a 
par,icu,ar bit refers to a subse. of points where .he va.ue of that bi. equals the 
desired value. The FALSE region is the compliment of the TRUE region. As an 
example, where i, is equal to zero, the TRUE region is the right-hand plane as 
shown in F.G. 3 and ,he FALSE region is .he left-hand plane. As another example, 
for the subse. of cons.e.,a.io„ poin, where bi, * equals zero the TRUE reg.on 
comprises ,he second and ,hird rows above and below the abscissa of ,he Cartes.an 
graph in FIG. 3 and ,he FALSE regions comprise the first and fourth rows above 
and below the abscissa. 

Addi,,onauy, the Karnaugh mapped QAM signals are also Gray coded such 
that ,he symbols adjoining a given symbo. only change by one bi, value. For 
example ,he symbo. 3W representing ,he bits 000000 is adjoined by four 
constellation poin, 302, 304, 306 and 308, two ,n the 1 direction and two » the Q 

directly below has the value 000001 where only bit q, has changed m value. 
The soft values of the demodulated symbols are formatted as a twos 
30 complement according to the following: 
I = sxxx.xxx = si./ 



where i refers to the integer portion bits, / is the fractional portion bits and s 
represents the sign bit. In mathematical terms I may be represented by the 
following expression: 

/= - .2"-' + I < ( 2' + f f,2-' - - s2»-> + I + / (8) 

A = 0 A: = 1 

for numbers represented by a N number of integer bits and a M number of 
fractional bits. The demodulated QAM symbol is also represented in the form of 
I+jQ- 

Particular properties and characteristics of the Karnaugh mapped QAM 
signal constellation illustrated in FIG. 3 afford simplified LLR calculations. One 
characteristic is that the constellation of FIG. 3 is mirror symmetric about the 
horizontal I and vertical Q axes for all bits except for the two most significant bits i, 
and q„ which respectively have opposing values on opposing halves of the Cartesian 
plane. It is noted that the symbol may be Karnough mapped differently such that 
other bits besides the two most significant bits have opposing values on opposing 
halves of the plane. However, typically only two bits are not mirror symmetric and 
these two bits are comprised of an i bit and a q bit. Another feature is that the log- 
likelihood ratio of a bit l } or is independent of the quadrature component of the 
demodulated symbol. Specifically, a bit \ } is independent of the Q quadrature 
component and a bit qj is independent of the I quadrature component. Thus, the 
calculation of the squared Euclidean distances as given in equation (7) previously 
can be greatly simplified since the calculation need only involve either horizontal or 
vertical distances, but not both. Hence, the log-likelihood ratio calculations for the 
Ij and Qj bits are given respectively by the following equations: 



LLR{ij)- min(A/]) 2 - min(A/°) 2 



(9) 



LLR(qj)= min(A2;) 2 - min(A0°) 2 (10) 



where AI S and AQj are the horizontal and vertical distances, respectively, 
between the demodulated symbol and a point where the bit equals a value that is 
either a zero or a one. 

Yet another property of the Karnaugh mapped signal that affords 
simplification of LLR calculations is that when a demodulated symbol resides in a 
FALSE region for a particular bit, the nearest point will lie in the row or column 
closest to the border between the TRUE and FALSE regions. Thus, the distance to 
the nearest point will be a value of one plus the magnitude of the distance from the 
demodulated symbol to the border. A value of one is added here since the borders 
in the Karnaugh map of FIG. 3 are set halfway between the rows and columns, 
which are spaced at a distance of 2. Hence, the distance from rows or columns to a 
boundary lying halfway between them is one. However, other embodiments may 
utilize distances between rows and columns other than two and, thus, different 
values would be added to the magnitude of the distance from the demodulated 
symbol to the border. Additionally, if the demodulated symbol lies in a TRUE 
region for a particular bit the distance to the nearest point is equal to fractional part 
/ or Independent on whether the integer portion i of the demodulated symbol 
component I is even or odd. 

Relying on these characteristics an algorithm is then formulated according to 
the teachings of the present invention whereby the number of calculations to 
compute the LLR of each bit is reduced. In the following discussion the disclosed 
algorithm is for a 64-QAM constellation. 

Using equation 9, discussed above, the log-likelihood ratio for the most 
significant bit /, is calculated using the expression: 



LLR{i x )= minCA/, 1 ) 2 - minCA/, 0 ) 2 



(10) 



Where I, is the twos complement formatted position of the first bit I,. In this 
case the minimum distance between the demodulated symbol bit I, and a 
constellation point having a bit I, with the value of zero can be determined using 



lh e renewing relationships rather Chan calculating scared Euclidean distances for 
all 32 symbols where the most significant bit will be equal to zero. 

First if the value of I is less than or equal to zero it is known that the 
modulated symbo, lies in the ,eft hand side of the Cartesian plane shown ,n HO . 
3 Thus it follows that the minimum distance to a constellation point havtng the 
bit equalto zero will be one (i.e., this value ,s arbitrary bu, is preferably one as 
illustrated in FIG. 3 where the first column of constellation po.nts are set at a value 
of one on the abscissa axis) plus the abso.ute value of I (i.e., the absolute value must 
be taken since I is less than or equal to zero or, in other words, negat.ve). 

Secondly, if the value of I is greater than zero the minimum distance to a 
constellation pom, having a bit ., equal to zero can be easily obtained dependent on 
whether the integer portion of the demodu,ated symbo, value . is even or odd. 
the integer portton is odd, then the minimum value is simply equal to the fracttona, 
portion of the twos complement formatted demodulated symbol. For example, tf 
the demodulated symbol location is at 3.9 in the i direction of the Cartestan plane 
then the Coses, constellation points will be those in the second column a, a value of 
3 Thus the distance along the I axis between the demodulated symbol and the 
constellation point is .9, which is the fractional part/of 1. Alternatively, if the 
integer portion , of the demodulated twos complemented symbol is even then the 
minimum va,ue is equal to 1-/. For exampie, if the demodulated symbo. lies 
in the i direction, the closes, constellation points will lie in the third courmn a, 5 
along the I axis of the Cartesian plane. Thus, the distance is 1 minus the fractional 
portion/(i.e., .l)or, in other words, .9. 

Finally if the demodulated symbo. lies beyond a value of 8 on the I axts of 
the Cartesian p.ane, then the minimum distance to a constel.ation symbo. w.th the 
most sign.ficant bit i, equa, to zero wil, be the difference between the value of 1 and 
the last column to the right of 7. However, since demodulated symbol values 
falling between 7 and 8 on the Cartesian plane are already accounted for by the 
previously discussed cond.tions, the distance is calculated by the expression 1 plus 

minus 8. . _ 

The foregoing conditions are set forth mathematically where - 0 as 
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follows: 



min(A/,°)= 1 + IM 
/ 

1-/ 
1+1-8 



I<0 

I > 0, Integer / is odd 
I>0, Integer f" is even 



Similarly, the minimum value of the distance between the demodulated 
symbol and the constellation point when the bit /, = 1 can be computed using the 
corollary conditions listed below. 



min(A/,')= 1 + |I| 
/ 

1-/ 
1+1-8 



I>0 

I<0, Integer i is odd 
I<0, Integer i is even 
K-8 



It is noted that for demodulated I symbol values falling on the Q axis, the 
algorithm arbitrarily is set such that zero is taken as a point in the right half plane. 
However, zero values could alternatively be assumed to be part of values in the left 
hand plane. 

Similarly, for finding the LLR of the bit i 2 , further characteristics of the 
Karnaugh mapped constellation illustrated in FIG. 3 are utilized to reduce the 
number of calculations necessary to derive the LLR. Here the i 2 bits are the same 
for the outer 2 columns on both sides of the plane whereas the inner 2 columns on 
each side of the Q axis are the same value for those particular constellation points. 
Accordingly, the cutoff points are the value of 4, which lies between the two 
regions on either side of the Q axis. Following are the specific equations used in 
the algorithm to determine the minimum distances for i 2 values equal to zero and 
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one, respectively, and also the equation used to determine the LLR of i 2 utilizing 
equation 9 above. 

Assuming i 2 = 0; 

min(A/ 2 °)= l+(|I|-4) |I|*0 

f |I| <4, Integer / is odd 

l_y |I| <4, Integer /' is even 

Assuming i 2 = U 

min(A/l)= l+(4-|I|) |I|<4 

/ |I|>4, Integer i is odd 

1 ./ J!) >4, Integer / is even 

l+(|I|-8) |I|>8 

Finally, to compute the last i bit i 3 , the particular characteristics of the 
Karnaugh mapped constellation are again utilized to reduce the calculation of 
squared Euclidean distances necessary to determine the LLR for this bit. Following 
are the conditions used to calculate the LLR for / 3 for the particular Karnaugh 
mapped 64-QAM constellation of FIG. 3 having the characteristics of the inner and 
outer columns on both sides of the Q axis having the same bit values for i 3 . 

Assuming i 3 = 0; 

min(A/ 3 °)= l+(2-|I|) |I|<2 

l+(|I|-6) |I|>6 

/ 2< 1 1|<6, Integer / is odd 

\.f 2< |I| <6, Integer / is even 



Assuming i : 



= 1; 



min(A/])= l + (|I|-2) 2<|I|<4 
l+(6-|I|) 4<|I|<6 
/ 

1-/ 
/ 

W 

l+(|I|-8) 



|I| <2, Integer / is odd 
|I| <2, Integer / is even 
6<|I| <8, Integer / is odd 
6<|I| <8, Integer i is even 
|I|>8 



Likewise, the Q, bits are also calculated using the same techniques that are 
used for the I, bits except that the distances are in the Q instead of I direction. 
Following is an example of the equations used in the algorithm for finding the LLR 
of the bit Q,. 



Assuming q l = 0: 

min(A0,°)= l + IQI 
/ 

1-/ 

l+Q-8 



Q<0 

Q>0, Integer q is odd 
Q>0, Integer q is even 
Q>8 



Assuming q x = 1: 

min(Ae, 1 )= 1 + 
/ 

1-/ 

l+Q-8 



Q>0 

Q<,0, Integer q is odd 
Q<0, Integer q is even 
Q<-8 



The present discussion omits discussing the calculation of the LLR's for bits 
Q 2 and Q 3 since they are corollaries to the calculations for bits i 2 and i 3 , 
respectively, discussed above. 

FIG. 4 illustrates a logic apparatus block diagram that may be used to 
5 implement the previously discussed algorithm. This apparatus may be employed 

within the log-likelihood ratio calculator 202 as illustrated in FIG. 2. In particular 
the demodulated symbol being formatted in the twos complement form of si./is 
input to the logic apparatus where the individual components are split into the sign 
s, the integer and fractional portions /./, the integer portion / alone and the 
10 fractional portion /alone. Prior to splitting the twos complemented symbol into its 

components, an absolute value block 402 is used to remove the sign portion S of the 
twos complemented symbol. Each of the values is input to a combinatorial logic 
404, which, in turn, calculates numerous flags dependent on the values input. For 
example, the logic 404 outputs a flag, when the value of I is greater than zero, an 
1 5 even or odd indication when the integer value is respectively even or odd or when 

the absolute value of I lies between, above or below particular numerical values. 
These flags are used by the logic to decide which of various inputs will be output, 
as will be discussed later and are essentially indicative of the position of I relative to 
various values along the I axis. 
20 The logic contains a number of summing blocks 406, 408, 410, 412, 414, 

416, 418, 420, 422, and 424 that all receive as an input the absolute value |I| of the 
received symbol. To this value, various values are either added or subtracted and 
the resultant values are used to input to the various decision multiplexers. Prior to 
input to the decision multiplexers, however, the resultants of the addition blocks are 
25 time controlled by D flip-flops being clocked with a clock cbl (i.e. , 426, 428, 430, 

432, 434, 436, 438, 440, 442 and 442). Hence, all of these values are computed in 
parallel simultaneously and are released to the decision multiplexers. Each of these 
decision multiplexers receives two or more input bits, which are comprised of the 
flag values output from the combinatorial logic 404. Furthermore, each of the 
30 multiplexers is used to output a value comprised of one of the inputs that represents 

the minimum distances AI J , AI J , AI \ , AI \ , AI \ , and AI \ . 



The decisions output by the multiplexers are based on the particular flag 
values input and correspond to the previously discussed conditions of finding the 
minimum distance based on the particular bit I, and whether the distances are to a 
zero or a one value. Of further note, the multiplexers are paired to output the AL 
values to both zero and one. For example, multiplexers 446 and 454 respectively 
output the values AI? and AI 1 , to D clocked flip-flops 448 and 456 that are 
triggered by clock cb2. As can be seen in FIG. 4, similar pairings of multiplexers 
(e.g., 462 and 470, and 478 and 486) are arranged, respectively, to determine the 
AI values for bits I 2 and I 3 )- The D flip-flops 448, 456, 464, 472, 480 and 488 are 
all simultaneously clocked to release the values of Af, to the zero or one values, 
respectively, to multipliers (i.e., 450, 458, 466, 474, 482 and 490) that are arranged 
to square the value that is input thereto. For each pairing of multipliers (i.e., 450 
and 458, 466 and 474, and 482 and 490) the squared resultants are output to 
respective summing blocks 452, 468 and 484 that subtract the squared minimum 
distances thereby calculating the LLR in accordance with the relationship described 
previously in equation 7. Each of the summing blocks 452, 468, and 484 are input 
to clocked D flip-flops, which are clocked in parallel via clock 4>3 to output the 
LLR values for the Ij bits simultaneously. It is noted that the apparatus used to 
calculate the LLRs of the Q component is identical with the apparatus illustrated in 



FIG. 4. 



An advantage of the configuration of FIG. 4 is that all of the operations are 
logically combinatorial thereby allowing execution of the operations to be clocked 
with a single clock and performed in a single clock cycle. Thus, the calculation of 
the flags output from the combinatorial logic 404 and the intermediate calculations 
output by the D flip-flops 426, 428, 430, 432, 434, 436, 438, 440, 442 and 444 
(and similarly those for the Q component) can be performed in a single clock cycle. 
Similarly, the minimum distances can be derived via the multiplexers and their 
accompanying flag inputs also on a single clock cycle. The multiplication operation 
performed by the multipliers configured to square the input value requires 2 clock 
cycles to perform. Finally the latched LLRs are released via the D flip-flops 460, 



476 and 492 in a single clock cycle from clock <$>3. 

Accordingly, the total number of clock cycles required to calculate the LLRs 
of a demodulated symbol is 5. Thus, the processing time for a 5376 symbol frame, 
assuming the same clock rate of 32 x 1.2288 MHz used previously, is 
approximately .684 milliseconds. The use of the above-described algorithm as 
implemented with the apparatus shown in FIG. 4, for example, reduces the time 
required to calculate the bit LLRs for a 5376 64-QAM symbol frame from 
approximately 150 milliseconds to less than 1 millisecond. Thus, feasible 
calculation of bit LLRs in 64-QAM modulated systems is possible even for high 
speed data channels. 

The above-described teachings of the present application achieve a quite 
significant reduction in the time required to calculate bit LLRs in a 64-QAM symbol 
communication system affording feasible use of bit decoding in high speed data 
communication channels. It is noted that the present method and apparatus do not 
apply to only 64-QAM but any M-ary QAM modulated system where the 
constellation of symbols may be square Karnaugh mapped. It will be further 
apparent to those skilled in the art, that other embodiments other than the specific 
disclosed embodiments described above may be devised without departing from the 
fair scope of the appended claims and their equivalents. 



